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ABSTRACT 

We present observations taken with the Donald C. Backer Precision Array for Probing 
the Epoch of Reionization (PAPER) of the Centaurus A field in the frequency range 120 to 
180 MHz. The resulting image has a dynamic range of 3000 and an rm.s. of 0.5 Jy beam^'. 
A spectral index map of Cen A is produced across the full band. The spectral index distri- 
bution is qualitatively consistent with electron reacceleration in regions of excess turbulence 
in the radio lobes, as previously identified morphologically. Hence, there appears to be an 
association of 'severe weather' in radio lobes with energy input into the relativistic electron 
population. We perform a detailed comparison of the large scale radio and X-ray emission 
from the ROSAT All Sky Survey. While the ROSAT field has significant gradients and struc- 
tures on 10° scales possibly unrelated to Cen A, two interesting correlations are seen between 
the radio and X-ray emission. First is an apparent 'cavity' generated by the northern radio 
lobe on a scale of 5°, possibly indicating excavation of thermal gas by the expanding radio 
source. Second is a correlation between radio and X-ray 'hot spots' at the end of the southern 
lobe, some 200 kpc from the nucleus. This likely arises from Inverse Compton scattering of 
the CMB by the relativistic electrons also responsible for the radio synchrotron emission. The 
magnetic fields derived from the (possible) IC and radio emission are of similar magnitude as 
the minimum pressure fields, ^ 1 /xG. 

Key words: galaxies: individual (Centaurus A) - instrumentation: interferometers - active: 
galaxies - galaxies: structure - X-rays: galaxies 



1 INTRODUCTION 



At a distance of 3.8 Mpc ( [Harris, Rejkuba, & Harris|20T0l l, Cen- 
taurus A is by far the closest radio galaxy, showing AGN-driven 
radio jet structures on scales from parsecs to 200 kpc jlsrael|1998] l. 
Given its proximity, and very large physical scale, Cen A has be- 
come a key tool in the study of radio source physics, and radio- 
mode feedback in suppressing late-time gas accretion in large el- 
liptical galaxies. However, its large angular scale (> 10°) presents 



a severe challenge to radio imaging of extended features with reas- 
onable spatial resolution. 

The last few years have seen a renaissance in low frequency 
radio interferometry, driven principally by the goal of studying the 
redshifted HI 21 cm emission from the neutral intergalactic me- 
dium during cosmic reionization over a redshift range of 6 < z 



< 14, or frequency range 180 MHz to 100 MHz (Fan, Carilli, & 
Keating|2006, FurlanetFo, Oh, & Briggs|2006; Morales & Wyithe 



2010 1. Numerous radio interferometers are currently in operation 
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to study reionization, including the Murcliison Widefield Array 
(MWA; Tingay et al. 20 12 ), the Low-Frequency Array (LOFAR; 
[Rottgering et al.||2006^ , and the Donald C. Backer Precision Ar- 
ray for Probing the Epoch of Reionization (PAPER; [Parsons et al.| 
|2010| >. In parallel to the reionization studies, the combination of 
good sensitivity, reasonable spatial resolution, and very wide field 
of view, of these interferometers presents a unique opportunity to 
image large-scale, non-thermal, radio continuum structures, such 
as those of Cen A. 

Here we present the imaging results on Cen A from the PA- 
PER array in South Africa. PAPER, with a field of view of approx- 
imately 60°, and a spatial resolution of --^ 20', is ideally suited to 
image Cen A. In section [2] we offer a brief overview of relevant 
studies of Cen A. In sections |3] and |4] we review the observational 
techniques and describe the routine used for reducing the PAPER 
data. We make a detailed comparison with the higher resolution 
morphological analysis of Cen A at 1.4 GHz by Feain et al. (201 1) 
in sectionjs] Section[6]presents the spectral index map of Cen A at 
148.2 MHz. We go on in section|7]to compare the radio and X-ray 
structures, and show that the spectral behaviour of Cen A is con- 
sistent with particle acceleration in turbulent structures in the radio 
lobes. Finally, section[8]summarizes our conclusions. 



2 CENTAURUS A 

Centaurus A is considered the archetype of Fanaroff-Riley type I 
radio galaxies (i.e. edge-darkened, lower luminosity radio galaxies 
- usually below 2 x 10^^ W Hz~ ^ sr~ ' , fpanaroff & Riley|| 1974| ). 
|Israel| ( (r998| l offers a comprehensive review of Cen A to that date. 

Cen A is part of the Centaurus A/M83 Group of galaxies. The 
host gala xy, NGC 5128, is a m assive elliptical galaxy (M = 4. 1 x 



van den Bergh|l976| , with a prominent, kpc-scale dust- 



10" Me 

lane. This dust lane likely indicates a recent merger with a dusty 
spiral galaxy, as is also suggested by the presence of optical shells 
and filaments around the host galaxy ( [Malin, Quinn, & Graham] 
|1983[[Schiminovich et al.| 1994 1. Stellar kinematics in the nuclear 
region of Cen A indicate a mass of the black hole of Mbh = (5.5 ± 
3.0) X 10^ Mq ICappellari etaTj ' 



2009 



At radio frequencies, early observations of Cen A revealed 
two giant outer lobes, extending up to 260 kpc north and south of 
the core, and possibly further (Coo per, Price, & Cole|1965| . Sub- 
sequent imaging has delineated radio structures from sub-parsec 
scale jets ( [Tingay et al.|1998y through two well defined inner lobes 
and jets on scales of ~5 kpc ( [Burns, Feigelson, & Schreier|1983[ ), 
as well as a 'middle lobe' in the northern part of the source on a 
scale of 30 kpc (Morganti et al. 1999V 

The most comprehensive radio imaging study of the outer 
lobes of Cen A to date was by [Feain et al.| ( (20TT| ). They combined 
observations by the Australia Telescope Compact Array and the 
Parkes 64 m radio telescope to image Cen A at 1.4 GHz with a 
spatial resolution of ~ 50" over a scale of 10°. In this image, they 
observe a number of new features in the outer lobes. The hook-like 
structure at the northern-most end of the northern lobe terminates in 
a series of radially-centric shells with decreasing brightness. These 
shells have similar curvature to a diffuse filament protruding from 
the eastern edge of the lobe, prompting [Feain et al.[j20lT| to sug- 
gest that both are modelled by the local group 'weather' {Safouris] 
[et al.|2009) . Near the northeastern edge of the hook they observe a 
partial ring and suggest that was created by the passage of a neigh- 
bouring galaxy. In the southern lobe, [Feain et al.[j2011^ identify a 
'vertex' and a 'vortex' close to 13*^ 20™ and —45° 15', which they 



show could be a region of plasma reawakened by strong shocks 
arising from the galactic core or by the passage of a dwarf galaxy. 
They also remark on the detection of a turbulent Faraday rota- 
tion measure signal jFeain et al.|2009) from the southern lobe and 
the presence of surface wave-like structures jBicknell, CameroTi^ 
|& Gingold 1990 ) close to the sharply defined western edge of the 
lobe and suggest that, together with the 'vertex' and a 'vortex', they 
could be the signature of Kelvin-Helmholtz instabilities. They ac- 
count for the 'vertex' and a 'vortex' being part of the wave-like 
structure by conjecturing a backflow dominated southern lobe with 
turbulent mixing. 

A detailed comparison of the X-ray and radio emission on 
scales of order 10 kpc shows clear evidence for the interaction 
between the expanding radio source and the thermal cluster gas 
( Ka rovska et al.|20 02 ). On much larger scales ('-^ 200 kpc), early X- 
ray imaging by the Ariel V satellite showed possible evidence for 
extended emission from the outer radio lobes ( |Cooke, LawrenceTj 
[& Perola|[T978t . [Cooke, Lawrence, & Perolaj (1978} suggest the 
emission is inverse-Compton upscattering of the CMB by the re- 
lativistic electrons in the radio lobes, and calculate a radio lobe 
magnetic field of 0.7 /iG, and that this field is comparable to the 
equipartition field. |Arp| ( [1994| ), presents the ROSAT sky survey X- 
ray emission on a scale of 10° around Cen A. This image shows 
a marginal excess of X-ray emission from the end of the northern 
radio lobe. Stawarz et al.[ ( [2012) present the latest X-ray analysis 
of Cen A with the Suzaku satellite. Their spectra of lobe regions 
show possible diffuse emission consistent with a thermal gas in the 
radio lobes with a temperature 0.5 keV and density ~ 10^^ cm^^, 
although the very large scales involved and the relative proximity 
of Cen A to the Galactic plane (situated 10° to 20° off the Galactic 
plane) leads to some confusion in the X-ray image. If real, this 
thermal gas has a pressure comparable to minimum energy. We re- 
turn to the X-ray-radio comparison below. 

Cen A has also been studied in the > 100 GeV energy range 
and its radio lobes are considered a likely origin for ultra-high en- 
ergy cosmic rays jFraija et al.|2012| ). 



3 OBSERVATIONS 

PAPER is an interferometric transit array aimed at detecting the 
fluctuations in the 21 cm emission from the Epoch of Reionization. 
The project currently consists of a 32-antenna test-array situated 
at NRAO, in Green Bank, USA, and a 64-antenna array in South 
Africa's Radio Quiet Zon^(latitude of approximately 30° S). The 
antennas are non-steerable, operate within a frequency range of 
100-200 MHz and have dual linear polarization capability. The 
primary beam is 60° FWHM with a spatially and spectrally smooth 
response near to the horizon ( [Pober et al.|2012) . 

We present here data taken during one of PAPER'S observing 
campaigns. The data were recorded in single polarization (E-W) on 
four consecutive nights (3"'/4* - 6*/7* of July 201 1) from approx- 
imately 4pm to Sam each 'night'. This observing campaign offers 
various investigation opportunities, such as studying the Galactic 
Centre or searching for radio transients and Supernova Remnants. 
Herein we present the imaging results on Cen A. 

At the time, the 64 antennas of the array were deployed in a 
minimum redundancy configuration ( [Parsons et al.|2012) within a 
300 m-wide circle, with an empty central circular region of radius 



' Further reference to PAPER will imply only the South African an'ay. 

© 0000 RAS, MNRAS 000, 000-000 



Imaging on PAPER: Centaurus A at 148 MHz 3 



150 



100 



50 



-50 



-100 



-150 



n r 




150 



100 



50 



-50 



-100 



-150 



Figure 1. PAPER «v-coverage for the full band and a 10-min dataset. The 
axes are labelled in units of wavelength. 



Table 1. Time intervals for Cen A observations. 



Day 


Time interval 




Number of 




(UT) 


(s) 


datasets 


030711 


16:12:35- 18:32:53 


8418 


15 


040711 


15:54:57- 18:25:16 


9019 


16 


050711 


16:09:09- 18:09:25 


7216 


13 


060711 


15:46:19-18:16:38 


9019 


16 



10 m. This generates a resolution of 26' at 145 MHz. Visibilities 
were integrated and recorded every 5.37 s. The full bandwidth is 
divided into 2048 channels, each of width 48.83 kHz. Since the 
antennas are not steerable, the observations were taken as a drift- 
scan of the sky. They were approximately continuous and divided 
into datasets of about 10-min long. Fig. [T| shows the Mv-coverage 
across the full PAPER band for one 10-min dataset. 

Only datasets in which Cen A was within approximately 1.5 
hours from zenith are considered, as over this hour angle, the 
primary beam gain is approximately constant along the source track 
(the primary beam's gain drops by about 32% between zenith and 
the outermost position). Table [T] shows which time intervals were 
used from each night and the number of 10-min datasets used per 
night. 



4 DATA REDUCTION 

Initial pre-calibration was performed using the Astronomical In- 
terferometry in PYthon (AIPY) packag^ This fixes hardware and 
header errors and subtracts cross-talk. An initial manual Radio Fre- 
quency Interference (RFI) excision step of the most prominently 
damaged channels was performed (mainly a few communication 
satellites bands, e.g. Orbcomm at 137 MHz). We find that due to 



paper's location, RFI only affects a few narrow channel ranges. 
Frequencies below 1 14 MHz and above 188 MHz were flagged due 
to bandpass roll-off. Within each 10-min dataset, we phase tracked 
the mean zenith, such that the result is similar to that of a tracking 
interferometer and can be imaged as such. 

The main calibration was carried out using the Astronomical 
Image Processing System (AIPS) packag^ A second round of 
manual RFI excision was performed by inspecting the scalar aver- 
ages for the total power spectra for each dataset and flagging outly- 
ing peaks. Approximately 67% of the initial data remained after 
both RFI excision steps. A final CLIPping at 3 times the average 
visibility level ensures that any remaining outliers are removed. 

Each 10-min dataset was divided into subsets of 102 frequency 
channels (4.98 MHz wide subsets) to mitigate bandwidth smearing 
effects. The 10-min time length of a dataset is sufficiently short for 
the effect of the source moving through PAPER'S primary beam to 
be negligible. All subsequent steps only used baselines longer than 
5 wavelengths, as the shortest baselines appear to be affected by 
very large scale structures which are not well sampled due to PA- 
PER'S non-uniform i(v-coverage (see Fig.[T|( at very short baselines 
(less than 10 m in length). 

Given the field of view of 60° of PAPER, the observed field 
cannot be approximated as a flat patch of sky anymore. We explored 
wide-field imaging in AIPS, which divides the field into a number 
of facets and cleans and images them quasi-independently. How- 
ever, we found that it does not produce a much improved image, 
but it does significantly increase the processing time. An altern- 
ative way of tackling this problem is by w-projection ( |CornwelUl 
|Golap, & Bhatnagar 2008 ), projecting the observed field to a tan- 



http://casper.berkeley.edu/astrobaki.index.php/AIPY 



gent plane, but this technique is not available in AIPS. 

The datasets were self-calibrated to obtain complex gain solu- 
tions for each antenna. Those with gain variations higher than 30% 
were flagged. Only four antennas had such high and unstable gains 
to warrant removal from at least one of the nightss. The following 
parameters were found to give good solutions after two phase self- 
calibration loops and a single phase and amplitude one: a solution 
interval of 2 minutes, a minimum number of 5 antennas allowed for 
a solution and a minimum allowed signal to noise ratio of 4. Clean 
components were used until the first negative. 

Each subset was CLEANed down to a set average flux level. 
This was taken to be the average root mean square noise (r.m.s.) of 
all the 10-min data subsets for a given night and channel rage and 
it varied by a factor of almost 10 between the lowest and highest 
channel range. Imaging was performed with a restoring beam of 
size of 1500" x 1300", the lowest resolution across the band. A uv- 
weighting scheme marginally naturaj^was found to be optimum in 
removing large scale effects, while still preserving enough inform- 
ation on the observed extended structures. 

The images were normalized by the primary beam model as 
a function of frequency as in the survey by [Jacobs et~aLl ( |20TT] l, 
and resampled onto the same sky-grid, using the datatset closest to 
zenith as reference. For each channel range, the images correspond- 
ing to different time intervals were averaged together. We measured 
the flux densities of bright sources within 18. 5"]^ from zenith for 
each channel range and manually identified them with Culgoora 
sources. We eliminated sources with inconsistent historical obser- 



' http://www.aips.nrao.edu/index.shtml 
* A ROBUST factor of - 1 in AIPS. 

^ This radius was chosen such that Cen A is included in the area thus 
defined. 
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Table 3. Peak flux densities for some sources in the Cen A field as found 
with the JMFIT task in AlPS from an average map of the middle thirteen 
channel ranges. The column names are self-explanatory; the expected flux 
is found from the Culgoora 160 MHz measurements and the corresponding 
spectral indices. 



Culgoora 


Expected 


Peak flux 


Rms 


source 


flux density 


density 


(Jy/beam) 


name 


(Jy) 


(Jy/beam) 




1221^23 


8.1 


7.2 


0.5 


1232^16 


11.0 


10.1 


0.5 


1252-289 


14.9 


17.1 


0.5 


1334-296 


2.8 


9.8 


0.5 


1355^16 


29.5 


31.9 


0.5 


1407^25 


9.6 


10.9 


0.5 


1424-418 


9.9 


8.3 


0.5 



vational data and obtained 13 sources to act as flux density calibra- 
tion. We fitted spectral energy distributions (SED) to the historical 
data for each of these sources with quadratic functions. The flux 
density calibrators, their positions and their derived SED paramet- 
ers are listed in Table |2] A least squares fit was performed for each 
channel range in order to derive scaling factors for the maps based 
on the calibrators' measured and expected SEDs. 



5 RESULTS 

Fifteen calibrated maps were produced, one for each channel range. 
Fig. [2] shows a weighted average of the middle thirteen. It has a 
dynamic range of approximately 3000 as found from the ratio of 
the r.m.s. of a region near the edge of the image and the peak flux 
of the core of Cen A. Its r.m.s. is 0.5 Jy/beam. Cen A dominates the 
centre of the image and the Galactic Plane crosses it further South. 
The Supernova Remnant (SNR) PKS 1209-52 is visible to the SW 
of Cen A, at —52° Dec, as a shell with two brightened limbs jDickel| 
|&Milne|I976l l. 

A lower limit of 5 wavelengths was imposed for calibration 
and imaging, which means structures larger than 1 1.5° are not im- 
aged. Very little data are eliminated this way and the negative dips 
in the maps evince the lack of short spacing data. Cen A is sur- 
rounded by a shallow region of negative flux. This is even more 
visible in and around the Galactic plane, since the Galaxy has ra- 
dio emitting structures on effectively all spatial scales. Again, these 
negative 'bowls' represent missing short-spacings in the interfero- 
metric images, and result in uncertainty in the total flux measure- 
ments of large sources, and diffuse surface brightnesses. Moreover, 
this effect gets worse with increasing frequency, and could bias the 
overall spectral index distribution analysis. Hence, the spectral in- 
dex discussion in section|6]will be mostly qualitative. 

However, there is some evidence that Galactic emission ex- 
tends all the way to Cen A (see the end of this section) and large- 
scale foreground emission will certainly have a significant galactic 
contribution. From |Landecker & Wielebinski| jI970^ , the bright- 
ness temperature in the direction of Cen A at 150 MHz is 360 K. 
Assuming an angular area of 12 deg^, the galactic foreground con- 
tribution over Cen A will be approximately 1 100 Jy. So, in a way, 
the 11.5° angular cutoff acts as a filter for part of this emission. 

Fig. [3] is a close-up view of Cen A. The core and inner lobes 
are saturated to show the structure in the outer lobes more clearly. 
The peak flux is found near 13^ 25™ 31' and -43° 01' 26" to be 
1631 ± 131 Jy. The integrated flux density over the entire source 
(taking the lowest contour level as the extent-limit for the lobes) is 



3792 Jy at the average frequency of 148.2 MHz of the map.'Haslam| 
|etal.|^I98I ) found the integrated flux density of Cen A at 408 MHz 
to be 2762 ± 370 Jy. Assuming an average spectral index over the 
entire source of -0.7 between 19.7 MHz and 960 MHz i jBoltonl 
& Clark||1960|l, then the expected flux at 148.2 MHz is 5611 Jy. 
Sheridan] (| I958[ l measured Cen A's integrated flux at 85.7 MHz to 



be 8700 ± 1300 Jy. Extrapolating again to 148.2 MHz gives 5929 
Jy. Both of these values are significantly higher than the one meas- 
ured here, stressing the importance of the missing short spaces, but 
neither [Haslam et"ar| t 1981 ) nor Sheridan ( 19581 correct their ob- 
servations for Galactic contributions, which, again, would amount 
to about 1 100 Jy. 

The northernmost part of the southern lobe displays emission 
from an extended background source, the radio galaxy MRC 1318- 
434B (Wright & Otrupcek 1990; Smith & Bicknell 1986), associ- 
ated with NGC 5090. This can be seen as a small bright patch on 
the northern edge of the southern lobe and dominates the meas- 
ured flux density. There is a second (compact) background radio 
source in the centre of the southern lobe, PKS 1320-446 l [Wright| 
|& Otrupcek|I990[ [Burgess & Hunstead||2006^ , but at our resolu- 
tion we estimate that this source is at most 30% of the flux density 
measured over a single synthesized beam, and does not affect the 
overall analysis below. 

The southern lobe shows a sharp edge along the southwestern 
boundary of the radio source, and more diffuse emission extending 
to the east in two 'prongs' of emission by about 2°. These two 
prongs are also seen at 1.4 GHz (Feain et al.). Also seen in this 
image is the 'gap' between the southern lobe and the nucleus, where 
the flux density drops below 7 Jy beam^ ' (as opposed to an average 
of 15 Jy beam^' on the northern edge of the southern lobe). This 
gap is also seen at 1.4 GHz. 

The northern middle lobe extends outwards from the inner 
lobe, parallel with it, and curves northward into the outer lobe, such 
that the outer lobe has a position angle of approximately 0°. At 
about 13** 28™ and —40° 33', the outer lobe curves east to form 
a hook-like structure. The northern lobe shows a general decline 
in flux density away from the core, with further brightening in the 
hook-like region at the northern end of the radio lobe, some 260 
kpc from the core. 

On the largest scales, faint emission may extended a few de- 
grees beyond the two-prongs to the southeast, although missing 
short spacings start to become an important issue in this regard. 
This very extended emission to the southeast has been noted pre- 
viously, and may relate to high-latitude Galactic emission, and not 
Cen A jBolton & Clark] 1 960| [Haslam et al.|I98"T] l. 



6 CEN A SPECTRAL INDEX MAP 

PAPER'S large bandwidth allows for spectral index maps of the 
observations to be computed across the band. The spectral indejj^ 
image is shown in Fig.|4] The spectral index was calculated at each 
point by taking the logarithmic difference between a weighted av- 
erage image of 4 lower frequency channels (123.58 MHz - 138.38 
MHz) and a weighted average of 4 higher frequency channels 
(158.11 MHz - 172.90 MHz). The r.m.s. of each image is used 
as a weight during averaging. Both images were blanked at a 5cT 
level. As mentioned before, incomplete reconstruction of the total 
flux will bias this calculation. 



* The convention 5 = v° is used. 
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Table 2. The 13 sources used as flux calibrators with their Culgoora name and positions and the derived SED fit parameters. The functional form logigS = 
A + B X logjQ V + C X logfg V. All coordinates are in J2000. 



Culgoora 
name 


RA 


Dec 




A 


B 


c 


Flux at 
148.2 MHz (Jy) 


1355-416 


13^ 55'" 57.4* 


-41° 38' 


14" 


8.0 


-0.8 





25.9 


1308-220 


13*' 08'" 57.6* 


-22° 00' 


33" 


-16.5 


5.0 


0.3 


67.3 


1407-425 


14^ 07" 23.5* 


-42° 32' 


56" 


15.9 


-2.6 


0.1 


10.9 


1413-364 


13'" 32.3* 


-36° 26' 


55" 


1.6 


0.5 


-0.1 


13.3 


1424-418 


14h 24m 45 9s 


-41° 52' 


58" 


16.9 


-3.3 


0.2 


8.2 


1421-382 


21" 11.6* 


-38° 13' 


29" 


9.0 


-1.0 





16.9 


1419-272 


14h 19111 54 5s 


-27° 14' 


13" 


22.2 


-4.0 


0.2 


21.9 


1422-297 


14h 22"! 32.9s 


-29° 47' 


06" 


13.5 


-2.1 


0.1 


17.0 


1206-337 


06'" 06.0* 


-33° 47' 


26" 


-3.2 


1.9 


-0.2 


12.8 


1218-373 


12^ 18'" 01.4* 


-37° 23' 


51" 


-8.2 


2.8 


-0.2 


8.0 


1221-423 


21" 04.0* 


-42° 19' 


30" 


-17.5 


4.6 


-0.3 


11.0 


1232-416 


12'' 32" 58.7* 


-41° 36' 


39" 


-9.6 


3.2 


-0.2 


9.7 


1252-289 


12'' 52" 00.7* 


-28° 57' 


13" 


5.6 





-0.1 


14.7 




1 10 



- a 



- -10 



Figure 2. PAPER final image of Centaurus A obtained with a 1500" x 1300" FWHM Gaussian restoring CLEAN beam. Sources listed in Tablejsjare marked 
with white circles. The two zoomed-in patches are on different greyscales to each other and to the general image. All coordinates are in J2000. 
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14'^00™ 13^^30™ 13'^00 



Figure 3. PAPER image of Centauras A obtained witli a 1 500" x 1 300" FWHM Gaussian restoring CLEAN beam. Contour levels are a geometric progression 
in \/2, staring from 2 Jy beam^' . Negative contours are in light grey. The green circles mark the positions of MRC 1318^348 andPKS 1320^46. Coordinates 
are in J2000. 



The northern lobe shows the spectral steepening over the first 
3° from the radio core, characteristic of FR I (edge-darkened) radio 
galaxies. Interestingly, the spectral index then becomes flatter in the 
northern hook-like structure, furthest from the core. 

The southern lobe also shows interesting spectral behaviour. 
The centre of the lobe and the southwestern hard-edge show flatter 
spectral indices, with steepening toward the more diffuse prong to 
the east. And again, the very end of the southern-most prong shows 
flattening of the spectrum towards its tip. 

The background radio galaxy on the northern edge of the 



southern lobe, MRC 1318-434B, can be seen as an unresolved flat- 
ter point in the spectral index map. 

7 ANALYSIS 
7.1 Spectral Index 

Spectral flattening in radio lobes is typically interpreted as a sign of 
energy input into the relativistic electron population ( |Carilli et al.| 
|1991[ l. While necessarily qualitative due to the short spacing issue. 
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Figure 4. Spectral index map of Centaurus A. Coordinates are in J2000. 



a comparison of the spectral index distribution and the morpholo- 
gical structures seen at high resolution, and their interpretation by 
|Feain et al.| ( |20n| , is illustrative. 

Starting with the southern lobe, we find that the central region, 
and the sharp southwestern edge, correspond to flatter spectral in- 
dices than the regions to the east (and close to the core, although 
this region is confused by a background radio source). This region 
is called-out explicitly in the Feai n et al.| ( |201 \\ analysis as be- 
ing rich in high contrast, turbulent structure, notably their 'vortex', 
'vertex', and 'waves'. 

Likewise, the end of the northern lobe was highlighted by 
[Feain et al.H201l] l as exhibiting multiple shells and a ring in the 
1.4 GHz image. This region also shows substantial spectral flatten- 
ing relative to the regions closer to the core. 

The spectral flattening at the southeastern edge of the southern 
lobe does not correspond to equally pronounced turbulent structure 
in the radio images, but it does overlap with some of the X-ray 
emission observed in the ROSAT image (see Section|T2|. 




10.30 



0.20 



_ - - 0.10 



Figure 5. X-ray image of Centaurus A, taken from the ROSAT all-sky sur- 
vey and smoothed with a 1500" x 1300" beam. Overlaid in white are the 
radio contours from IJy up to 100 Jy. Coordinates are in J2000. 



7.2 ROSAT All Sky Survey Comparison 

Fig. |5] shows a greyscale of the 0.1 - 2.4 keV emission from the 
ROSAT All Sky Surve3Q(RASS) with the PAPER Cen A contour 
overlaid. The Xray image is admittedly noisy, with large gradients 
and clumpy structure, however, there are two interesting correla- 
tions between the radio and X-ray to note, beyond the well studied 
correlation on 10 kpc scales of the radio and X-ray emission from 
the core and inner jets. 

First, in the northern lobe there appears to be a deficit of X-ray 
emission due to the radio lobes. This is most clearly seen along the 
eastern edge of the lobe. Such X-ray deficits in clusters and group 
X-ray halos are often seen associated with cluster radio galaxies, 
and likely correspond to an evacuation of the thermal cluster gas 
via the expanding radio source (McNamara & Nulsen 2007; ,Fabian| 
|2012[ l. This cavity would then provide a relative pressure probe of 
the cluster gas, given the synchrotron and IC constraints below. 
Unfortunately, the large X-ray filament running through the field 
confuses this signature. It is also possible that this large scale fea- 
ture (scale of tens of degrees), running roughly north-south could 
be associated with high latitude Galactic emission (Stawarz et al.| 
|2012[ l, leading to a false deficit detection. This needs to be con- 



www.xray.mpe.mpg.de/cgi-bin/rosat/rosat-survey 
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13''35' 



Figure 6. Close-up of the southern lobe: radio contours in white overlaid 
on a high pass filtered X-ray greyscale with the mean surface brightnesses 
at 4 times lower spatial resolution subtracted to show high contrast features. 
Coordinates are in J2000. 



Table 4. Magnetic field strengths found from the RASS X-ray image. 
Lx-ray is the measure X-ray luminosity, Bi^jnp the magnetic field strength 
found through the minimum pressure condition and Bjc the magnetic field 
strength found through assuming inverse Compton scattering. S 1 and S2 are 
marked in Fig.[6]and a is the spectral index. We use two values for the spec- 
tral index to account for the uncertainty in the absolute a measured from 
Fig.g 







SI 




S2 




a = -l 


a = -1.8 


a = -l 


a = -1.8 


log(Lx-ray(erg s ')) 


39.0 


39.2 


39.9 


40.1 




0.56 


0.80 


0.67 


0.96 


Bic (MG) 


0.25 


1.25 


0.12 


0.74 



between 0.1 jlG and 1.3 /iG, again depending on the spectral index. 
In general, the fields based on a minimum pressure assumption and 
those based on inverse-Compton emission are of similar magnitude. 

[Hardcastle eTar] p009| l find an upper limit for the density 
of the X-ray emitting gas in Cen A's lobes of lO^'* cm^^, which 
agrees with the value range derived by |Stawarz et al.| ( |2012[ l. As- 
suming Ug = 10^^ cm^^ and a gas temperature of kT ~ 0.5 keV (as 
found by |Stawarz et al.|2012^ , then the lobes' thermal pressure is 
approximately 8 x 10^ erg cm^ . Stawarz et al. ||2012l| found that 
this thermal pressure counterbalances the pressure of a non-thermal 
plasma (8 x 10^'^ erg cm^-*) with an average magnetic field of 
strength B --^ 1 in rough equipartition with ultra-relativistic 
electrons. It is encouraging that the magnetic field strength required 
for pressure balance of the different plasma species is of the same 
order as the average magnetic field strength we find for the Cen A 
lobes, B ~ 1 ;tiG. 



firmed through future large scale X-ray imaging of Cen A. If the 
deficit is indeed a cavity due to the lobe of Cen A, this would be 
the first direct evidence for large scale X-ray emission from the Cen 
A group. 

More noticeable is the correlation between a series of X-ray 
and radio knots in the southern end of the radio source. Fig. |6] 
shows a close up of the region. Two X-ray knots are seen, follow- 
ing closely the curvature of the southern lobe and radio knots to 
the southeast. We have performed a pixel-pixel correlation in this 
region and find that the peaks correlate reasonably well, especially 
SI (see Fig.|6j correlation coefficient of 0.94). Again, the ROSAT 
image is somewhat confused, but the correlation of the peaks and 
curvature of the features suggest real association. 

We calculate the magnetic field strength in these regions in 
two ways. Fist, we consider minimum pressure conditions ( |Miley| 
|1980j ), with a relativistic electron to proton ratio of k = 1, unity 
filling factor, upper and lower frequency cutoffs for the radio spec- 
trum of 10 MHz and 10 GHz and a luminosity distance of 3.75 
Mpc. We assume that the emission is coming from a spherical re- 
gion in the lobe of 15' radius. The resulting fields are given in Table 
|4] The values range from 0.5 /xG to 1 /iG, depending on the spectral 
index assumed. 

Second, we consider that the X-ray emission is due to inverse 
Compton upscattering of the ambient CMB photons (the dominant 
photon field in the outer lobe regions) by the synchrotron emitting 
relativistic electrons jHarris & Grindlay||T979l |Miley|[T980) . The 
X-ray luminosity was determined after convolving the counts to 
15' resolution, and subtracting a background based on an annulus 
with the same area as the photometric circle. The resulting X-ray 
luminosities and magnetic fields are given in Table|4] We find fields 



8 CONCLUSIONS 

We present wide field imaging of large scale structures taken with 
the Precision Array for Probing the Epoch of Reionization in the 
frequency range 1 20 to 1 80 MHz of the closest radio galaxy, Cen- 
taurus A, and a comparison to the large scale X-rays. 

We observe a correspondence between regions rich in com- 
plex morphological radio structure at 1.4 GHz (described by |Feain| 
|etal.|2011^ and regions of flatter spectral index in our 148.2 MHz 
maps. The northern loop of the northern lobe, where |Feain et al.] 
(2011) show a series of shells and a ring, sees a flattening of the 
spectral index relative to the regions closer to the core. Similar flat- 
tening is seen in the central and western part of the southern lobe, 
where [Feain et al.] ( |2009^ observe a 'vertex' and a 'vortex'. This 
likely indicates active particle acceleration on large scales in the 
radio lobes, through supersonic turbulences. We note that the gap 
between the core of the galaxy and the northern edge of the south- 
ern lobe, which has also been observed at 1.4 GHz jFeain et al.| 
[20TT] l, is still seen at 148.2 MHz. 

Considering the RASS comparison, the end of the southern 
lobe shows a correlation of X-ray emission with a series of radio 
knots seen in the 148.2 MHz PAPER image. Assuming this is in- 
verse Compton up-scattering of the CMB by the synchrotron emit- 
ting relativistic electrons, we find that the magnetic field strengths 
required are comparable to the ones found from the ROSAT image, 
assuming a spherical emission region embedded in the southern 
lobe. Also, fields of the order of 1 jlG are found across both lobes 
from the radio map assuming minimum energy conditions. Assum- 
ing an average minimum energy of ultra-relativistic electrons of 
4 X 10^'^ erg cm^-* and cylindrical symmetry for Cen A (radius 2° 
and height 4°), we find the total non-thermal energy stored in the 
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lobes to be 1.7 X 10 erg. The possible cavity seen in the northern 
lobe due to the expanding radio source provides a rough measure 
of the thermal pressure in the cluster, via balance with minimum 
energy pressures in the lobes. 
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